e+e 3 jets and event shapes at NNLO 
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We report on the calculation of NNLO corrections to the 3-jet cross section and related event shape distributions 
in electron-positron annihilation. The corrections are sizable for all variables, however the magnitude of the 
corrections is substantially different for different observables. We observe that inclusion of the NNLO corrections 
yields a considerably better agreement between theory and experimental data both in shape and normalization of 
the event shape distributions in the region where the perturbative result is expected to hold. A new extraction of 
Oa using the event shape variables up to NNLO yields a considerably better consistency between the observables 
indicating a stabilization of the perturbative corrections at this order. 



1. Introduction 

Jet observables in electron-positron annihila- 
tion play a pivotal role in studying the dynam- 
ics of the strong interactions, described by the 
theory of quantum chromodynamics (QCD). In 
addition to measuring multi-jet production rates, 
more specific information about the topology of 
the events can be extracted using variables which 
characterize the hadronic structure of an event. 
With the precision data from LEP and SLC, 
experimental distributions for such event shape 
variables have been extensively studied [112] and 
have been compared with theoretical calculations 
based on next-to- leading order (NLO) parton- 
level event generator programs [3|4j , improved by 
resumming kinematically-dominant leading and 
next-to-leading logarithms (NLO-t-NLL) j5] and 
by the inclusion of non-perturbative models of 
power-suppressed hadronisation effects [6] . 

Up to now, the precision of the strong cou- 
pling constant determined from event shape data 
has been limited largely by the scale uncertainty 
of the perturbative NLO calculation. We report 
here on the first calculation of NNLO corrections 
to the 3-jet cross section and related event shape 



variables. The knowledge of the NNLO correc- 
tions to the event shape distributions has impor- 
tant phenomenological impact on the extraction 
of as from LEP data. 

2. The 3-jet cross section at NNLO 

Jets are defined using a jet algorithm, which 
describes how to recombine the momenta of all 
energetic hadrons in an event to form the jets. 
These algorithms are used in the experimental 
analysis and in the parton-level event generators 
to combine particles into jets. Here we present 
the first calculation of the NNLO corrections to 
the three-jet production rate at parton-level in 
e+e^ annihilation using the Durham measure [7]. 

The calculation of the corrections for three- 
jet production is carried out using a recently 
developed parton-level event generator program 
EERAD3 [8J which contains the relevant matrix el- 
ements with up to five external partons. Besides 
explicit infrared divergences from the loop inte- 
grals, the four-parton and five-parton contribu- 
tions yield infrared divergent contributions if one 
or two of the final state partons become collinear 
or soft. In order to extract these infrared diver- 
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Figure 1. Perturbative fixed-order description of 
the three-jet rate in the Durham jet-scheme at 
Q = Mz, compared to data obtained with the 
ALEPH experiment j2j. 



gences and combine them with the virtual cor- 
rections, the antenna subtraction method 9^ was 
extended to NNLO level [10] and implemented for 
e+e" 3jets [TT and related event-shape vari- 
ables [n] into EERAD3. 

Figure [1] displays the three-jet rate at LEPl en- 
ergy Q — Mz as function of the jet resolution j/cwt 
at LO, NLO, NNLO. The theoretical uncertainty 
band is defined by varying the renormalization 
scale fj, in the coupling constant in the interval 
Mzf^ < < 2Mz, and the average value [13] 
as{Mz) — 0.1189 is used, consistently evolved to 
other scales at each order. 

Since the error band in the region 10^^ > 
Ucut > 10~^ is barely visible in the plot, we dis- 
play the relative theoretical uncertainty 

g _ max^(cr(^)) - min^(fT(/i)) 
2a{^l = Mz) 

at NLO and NNLO as an inset. The uncertainty 
on the LO calculation is constant at 10.2%. 

As can be seen from the plot, the theoretical 
uncertainty is lowered considerably compared to 
NLO. Especially in the region 10~^ > Ucut > 
10^^, which is relevant for precision phenomenol- 
ogy, one observes a reduction by almost a factor 
three, down to below two per cent relative uncer- 



tainty. 

For large values of Ucut, Ucut > lO^'^, the NNLO 
corrections turn out to be very small, while they 
become substantial for medium and low values 
of Ucut- The maximum of the jet rate is shifted 
towards higher values of ycut compared to NLO, 
and is in better agreement with the experimental 
observation. 

The fixed-order theoretical predictions for the 
three-jet rate become negative for small values 
of 

Ucut, where fixed order perturbation theory is 
not applicable due to the emergence of large loga- 
rithmic corrections at all orders which require re- 
summation |7|14j . We therefore restrict our com- 
parison to Tjcut > 10~^. Even with this restric- 
tion, at low jet resolution, the fixed-order NNLO 
description lies above the data. The theoretical 
parton-level prediction is compared however to 
hadron-level data, thereby neglecting hadronisa- 
tion corrections, which may account for part of 
the discrepancy. 

The total hadronic cross section consists of the 
sum over all jet multiplicities. At ©(a^), this sum 
runs from two-jet through to five-jet final states, 
such that the corresponding fractional jet rates 
must add to unity. Consequently, our calculation 
yields the N'^LO expression for e^e~ — > 2 jets as 
a by-product. 

Figure [2] shows the parton-level theoretical pre- 
dictions for the jet fractions at first, second and 
third order in the strong coupling constant, com- 
pared to experimental hadron-level data from 
ALEPH [2]. 

By comparing the three plots, we observe that 
the agreement for each of the jet rates becomes 
systematically better as the order of perturbation 
theory increases. At each order a new multi-jet 
channel opens up, e.g. the five-jet rate at 0{al), 
which is positive definite and essentially mono- 
tonically increasing as ycut decreases. Since all 
jet rates are normalized to unity, the new five- 
jet channel has the effect of reducing the con- 
tribution to the two-jet, three-jet and four-jet 
rates, in the region of log]^o(2/ciit) where the five- 
jet rate contributes. One very clear effect is to 
cause the turnover in the four-jet rate (which is 
not present at 0{a1)). A second effect is to add 
more structure to the shape of the two- and three- 
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Figure 2. Jet rates at first, second and third order in the strong couphng constant, compared to data 
from ALEPH [2]. The rates are normahzed to the total hadronic cross section at that order. 



jet rates, which he much closer to the data for 
logio(ycMt) < —2.5. Of course, the effect of the 
higher order corrections also extends to larger val- 
ues of ijcut and, by adding more structure to the 
theoretical prediction, one obtains a better de- 
scription of the data. 

3. Event shape variables 

In order to characterize hadronic final states in 
electron-positron annihilation, a variety of event 
shape variables have been proposed in the litera- 
ture. For a review see e.g. [15 16 . In our study 
we considered only variables for three-particle fi- 
nal states which are thus closely related to three- 
jet final states. Among these event shapes, six 
variables were studied in great detail [T]: the 
thrust T, the normalized heavy jet mass p, the 
wide and total jet broadenings Bw and Bt, the 
C-parameter and the transition from three-jet to 
two-jet final states in the Durham jet algorithm 
called Y3. 

The perturbative expansion for the distribu- 



tion of a generic observable y up to NNLO at 
e~^e~ centre-of-mass energy ^/s, for a renormal- 
ization scale fi^ involves perturbative coefficients 
|12j which only depend on the event shape vari- 
able y itself and the strong coupling constant 
as- Those coefficients are computed by a fixed- 
order parton-level calculation, which includes fi- 
nal states with three partons at LO, up to four 
partons at NLO and up to five partons at NNLO. 

The precise size and shape of the NNLO correc- 
tions depend on the observable in question. Com- 
mon to all observables is the divergent behaviour 
of the fixed-order prediction in the two-jet limit 
(corresponding to small values of the event shape 
variable y), where soft-gluon effects lead to an en- 
hancement of the fixed-order coefficients by pow- 
ers of ln(l/?;). In order to obtain reliable predic- 
tions in the region of y <C 1 it is necessary to 
resum entire sets of logarithmic terms at all or- 
ders in as- A detailed description of the predic- 
tions at next-to-leading-logarithmic approxima- 
tion (NLLA) can be found in Ref. [I6] . 

For several event shape variables (especially T 
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Figure 3. Thrust distribution at Q = Mz at LO 
(blue), NLO (green) and NNLO (red). The sohd 
hnes represent the prediction for renormahsation 
scale ji = Q and as{Mz) — 0.1189, while the 
shaded region shows the variation due to varying 
the renormahsation scale between jj, = Q/2 and 
jjL = 2Q. The data is taken from [5]. 



and C) the full kinematical range is not yet re- 
alised for three partons, but attained only in the 
multi-jet limit. For the thrust distribution 1 — T, 
as seen in Fig. [3l the multi-jet limit corresponds 
to 1 — T > 0.5. Consequently, the fixed-order de- 
scription is expected to be reliable in a restricted 
interval bounded by the two-jet limit on one side 
and the multi-jet limit on the other side. 

In the intermediate region, we observe that in- 
clusion of NNLO corrections (evaluated at the Z- 
boson mass, and for fixed value of the strong cou- 
pling constant) typically increase the previously 
available NLO prediction. The magnitude of 
this increase differs considerably between differ- 
ent observablesfT^. it is substantial for T (18%), 
Bt (17%) and C (15%), moderate for p and Bw 
(both 10%) and small for ^3 (6%). For all shape 
variables, we observe that the renormalization 
scale uncertainty of the NNLO prediction is re- 
duced by a factor 2 or more compared to the NLO 
prediction. We observe that the NNLO prediction 
describes the shape of the measured event shape 
distributions over a wider kinematical range than 
the NLO prediction, both towards the two-jet and 



the multi-jet limit. 

4. Determination of the strong coupling 
constant 

Using the newly computed NNLO corrections 
to event shape variables, we performed [17] a 
new extraction of as from data on the standard 
set of six event shape variables measured by the 
ALEPH collaboration [2^ at centre-of-mass ener- 
gies of 91.2, 133, 161, 172, 183, 189, 200 and 206 
GeV. The combination of all NNLO determina- 
tions from all shape variables yields 

as{Mz) = 0.1240 ± 0.0008 (stat) 
± 0.0010 (exp) ± 0.0011 (had) 
± 0.0029 (theo). 

We observe a clear improvement in the fit qual- 
ity when going to NNLO accuracy. Compared 
to NLO, the value of as is lowered by about 
10%, but still higher than for NLO+NLLA 0, 
which shows the obvious need for a matching of 
NNLO-I-NLLA for an even more precise result. 
Work is in progress in this direction [18] . 

As can be seen in Figure [H the scatter among 
the ck^-values extracted from different shape vari- 
ables is lowered considerably when going to 
NNLO accuracy and the theoretical uncertainty 
is decreased by a factor 2 (1.3) compared to NLO 
(NLO-fNNLA). The different sizes of the NNLO 
corrections for different observables is responsi- 
ble for these large improvements. One infers that 
the scatter present at next-to-leading order was 
largely due to missing higher order perturbative 
corrections. 

5. Outlook 

Our results for the NNLO corrections to the 
3-jet cross section and related event shape distri- 
butions in electron-positron annihilation open up 
a whole new range of possible comparisons with 
the LEP data. The potential of these studies is 
illustrated by comparisons of the NNLO fixed or- 
der results with jet and event-shape data from 
ALEPH. The corrections are sizable for all vari- 
ables, but yield a considerably better consistency 
between the observables indicating a stabilization 
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Figure 4. The measurements of the strong coupling 
constant as for the six event shapes, at ^/s = M^, 
when using QCD predictions at diflerent approxima- 
tions in perturbation theory. The blue band indicates 
the uncertainty due to renormalisation scale variation 
in each theoretical description. 



of the perturbative corrections at this order. A 
fit to event shape data yielded a new determina- 
tion of as. We anticipate a further improvement 
by matching of the fixed order NNLO calculation 
with NLLA resummations jl8j . 
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